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ABSTRACT: We compute the zero frequency current noise numerically and in several 
limits analytically for the coulomb blockade problem consisting of two tunnel junctions 
connected in series. At low temperatures over a wide range of voltages, capacitances, and 
resistances it is shown that the noise measures the variance in the number of electrons in 
the region between the two tunnel junctions. The average current, on the other hand, only 
measures the mean number of electrons. Thus, the noise provides additional information 
about transport in these devices which is not available from measuring the current alone. 
1992 PACS 72.70.+m , 73.40.Gk 



I. INTRODUCTION 



When the charging energy of a tunnel junction is larger than the temperature, elec- 
tron tunneling events across the junction become correlated. These correlations lead to a 
variety of phenomena which fall under the rubric of single electron charging effects and 
the Coulomb blockade. ''' Recently, there has been renewed interest in the Coulomb 
blockade because of the wealth of new physical realizations: metal-insulator-metal tunnel 
junctions with small metal particles in the insulator, ^'^'^ lithographically patterned tun- 
nel junctions,^ scanning tunneling microscopy of small metal droplets,^ narrow insulating 
wires, -"^^ and even thin crossed wires. ■'^■'^ In some cases the technology has advanced to such 

an extent that one can consider making practical devices based on the coulomb blockade. 
12,4 

The theoretical work in this area has focused on the average current in either DC or 
AC measurements. There has been some work examining the effect of noise in the external 
circuit on the average current;^^ however, there have not been any studies until recently 
treating the current noise as an interesting phenomena in itself. ■'^'^'^^'-'^^ In this paper we 
compute the zero frequency current noise in one of the simplest coulomb blockade devices 
consisting of two tunnel junctions connected in series. We apply the same master equation 
used to compute the average current.'^'-'^^'-'^^ This work is a direct outgrowth of an earlier 
paper in which we computed the noise with a similar equation, which did not have charging 
effects explicitly. ■'■^ 

We study the zero frequency as opposed to the finite frequency noise because the time 
scale at which the frequency dependent noise begins to show structure is the "RC" charging 
time of the system, which is quite small (~ 10~^ sec).^) Noise experiments, on the other 
hand, are typically done in the regime below 10^ sec~^. Also, because this calculation 
is based on the usual master equation for the coulomb blockade, it will contain thermal 
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noise and shot noise, but not 1/f noise, which is often seen in tunnel junctions at low 
frequencies.^^ 1/f noise in tunnel junctions is usually assumed to be due to defects in or 
near the junction. These defects are not contained in our model. Although 1/f noise will 
dominate at the lowest frequencies, there is often an intermediate regime between the very- 
low frequency regime and the high frequency regime where shot noise dominates. By 
focusing on the thermal and shot noise as opposed to the 1/f noise we are able to make 
quantitative predictions for the noise. 

The zero frequency current noise is interesting for several reasons. First, it is a measur- 
able effect. It may limit the accuracy of some average current measurements. Secondly, we 
shall see that for this problem the noise provides a measure of the variance in the charge 
fluctuations as a function of time. The average current only provides information about 
the average charge between the two capacitors. Finally, the noise versus current curve 
has a rich structure, which is particularly revealing when rescaled by the average current. 
This can be used either to fit the parameters in the model more accurately or to provide 
a consistency check for parameters already determined from the average current. Thus, 
noise measurements provide an additional test for the underlying rate equation used to 
describe transport in this simple coulomb blockade device. 

This paper is organized into three parts: formalism (Sec. II.), analytic results (Sec. 
III.), and numerical results (Sec. IV.). In the formalism section the model and underlying 
master equation are reviewed, the formulas for the average current and noise are intro- 
duced, and the numerical technique used to compute the current and noise is discussed. 
Next in the analytic results section we calculate the noise in four limits: (i) the zero voltage 
limit, (ii) the region where the voltage is larger than the temperature but still smaller than 
the charging energy, (iii) the region just after the onset of current flow, and (iv) the large 
voltage limit. The details of the fluctuation dissipation theorem (zero voltage limit) for 
this system are presented in the appendix. Finally, the noise is calculated numerically and 
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compared to the analytic results. All the results are summarized in Sec. V. 



II. FORMALISM 



A. System 

The system which we study consists of two tunnel junctions connected in series (see 
Fig. 1(a)). The junctions are denoted by L for left and R for right, and have resistances 
Rl and Rji and capacitances and Cji, respectively. We also introduce the net resistance 
and capacitance of the junctions: R = R^ + Rr and C = + Cr. The capacitance C 
is the capacitance which enters the charging energy of the junctions, Eq — e^/(2C), even 
though the capacitance of the two capacitors in series is {C^^ + C^^)~^. To drive the 
current a voltage V is applied to the left junction, while the voltage on the right junction 
is set to zero. The voltage in the central or middle region between the two junctions, Vm-, 
fluctuates depending on the number of excess electrons in this region, n. The voltage drops 
across the left junction, — ^^(n), and the right junction, Vm(^) — are found using 
classical electrostatics to be 

V,-VM{n) = ^V+'^ + V, (1) 
VM{n)-Vn^^V-'^-V,. (2) 

Here and in all following discussions the magnitude of an electron's charge is e so that 
— e is the electron's charge. The additional voltage Vp has been included to account for 
any misalignment of the Fermi level in the middle region with the Fermi levels of the left 
and right leads when V and n are zero. ^'^^ One can also include an external gate on the 
central region, ■'^^ which has a similar effect to Vp. For simplicity we omit such a gate. 

4 



B. Master Equation 



The transport through the tunnel junctions is governed by four tunnehng rates: the rate 
for electrons to tunnel onto the central region from the left (T^^^^i) and right (r^^^_|_-|^) 
and the rate for electrons to tunnel off of the central region to the left {T^^^_^) and right 
(r^^„_l).^'^'^'^^ The number of excess electrons in the central region between the two 
junctions is n. These rates are computed via Fermi's golden rule. In order to write the 
rates we introduce a function 7 (e): 

For energies e larger than the temperature, hgT, 7 (e) is approximately e, showing that 
the tunneling rate increases as the allowed phase space is increased. For energies e less 
than —ksT, 7(e) is exponentially suppressed: 7(e) ~ lele"^!*"', since the energy for these 
processes must come from thermal fluctuations. In terms of the charging energy, Ec = 
e^/(2C), the tunneling rates are^'-^^'-^^ 

TnSnil = 7 (±e[lV(n) - V,^^,] - Ec) . (4) 

The energies which enter the rates of Eq. (4) are the voltage drops offset by the charging 
energy, indicating that tunneling is suppressed for voltages smaller than the charging 
energy. 

The state of the system at time t is described by the probability pn{t) that there are n 
electrons in the central region. Clearly the sum of the pn{t) is unity: J2nPn{t) ~ ^- '^^^ 
time evolution of pn{t) is governed by the net rate, F^^j = ^f_^j + Tfl^j, to go from i to 
j excess electrons in the middle region. By incorporating the rates Fj^j in a matrix M, 

1 0, otherwise, 
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the Master equation for the time evolution of p{t) may be written succinctly as 

A direct consequence of Eqs. (5) and (6) is that the sum of the Pn{t) is independent of 
time because 

J^Mij^O. (7) 

i 

Equation (7) in turn implies that the matrix M has a zero eigenvalue so there is a steady 
state solution, i.e., a vector p^°^ which satisfies 

= Mp(°). (8) 
C. Current 

Although only the net transition rate, Tj^j, enters in computing p{t), the current 
depends on whether an electron travels to the right or left. For example, a transition 
from n to n + 1 electrons in the central region has the rate T^^n+l- If the additional 
electron tunneled from the right, this process gives a negative contribution to the number 
current and a positive contribution to the electrical current. On the other hand, if the 
additional electron tunneled from the left, there would be a positive contribution to the 
number current and a negative contribution to the electrical current. To take into account 
the difference in sign for the two processes, we introduce two new matrices, and v^, 
which contain the rates for tunneling across the left and right junctions, respectively. The 
sign of a matrix element is positive if a process gives a positive contribution to the number 
current from left to right and negative if a process gives a negative contribution to the 
number current. 



V- • 



if^ = i-l, 

H-)r',Z ifi = j + i (9) 

0, otherwise. 
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With these definitions the electrical current across the left (II) and right (Ir) junctions at 
time t are 

hiuM = -e^(v^(^V(t)), ^ -eTVK(«)p(t)}. (10) 

i 

In Eq. (10) we have defined a trace of a vector to be the sum of its elements. Both the 
currents 1^ and Iji are positive when the electrical current goes from left to right. To 
obtain the steady state current from Eq. (10) p{t) is replaced by 

In general the current in the left and right junctions will be different; however, in the 
zero frequency limit they must be the same because of the continuity equation. To derive 
the continuity equation we introduce the number matrix N, where N^j — i x S^j. The 
expectation value for the number of electrons in the central region at time t is 

N{t) = TV{Np(t)}. (11) 
Using the rate equation for the probability, Eq. (6), the time derivative of N{t) is 

^ = TV{NMp(t)}. (12) 
The key observation is that — is the commutator of N and M. 

v^-v^=[N,M]. (13) 

Using Eqs. (7) and (1.14) the time derivative of the charge in the central region may be 
written as 

= -eTr{[N, M]p{t)} = h{t) - /^(t), (14) 

which is the continuity equation. As noted above in steady state, dN{t)/dt = 0, and II is 
equal to Iji. 



D. Noise 



To define the noise we introduce the propagator P{t) which gives the time evolution 
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of pn{t) 

P(t) = exp(Mt). (15) 

The conditional probability that there are m electrons in the middle region given that 
there were n electrons at t = is Pm,nit)- It is understood that t is positive in Eq. (15). 
With P(t) we can compute all possible correlation functions between the density and the 
current. ^^'^"^ For example, the density-density correlation function is 

{N{t)N{Q)) = ^(t)Tr{NP(t)Np(°)} 

(16) 

+ ^(-t)T>{NP(-t)Np('')}. 
This equation has a simple physical interpretation. Initially, the probability distribution 
is the steady state solution, p^°\ The number of the electrons in the central region is 
measured by N at this initial time. Next, the system is propagated forward in time via P, 
and the number is again measured with N. For t > the initial time is t = and the final 
time is t, while for t < the initial time is t and the final time is 0. In a similar manner 
the density-current correlation function is (^ = i? or L) 

{N{t)lA{0)) = - e^(t)Tr{NP(t)vV°^} 

(17) 

-e^(-t)Tr{v^P(-t)Np(°)}. 
The only difference between Eqs. (17) and (16) is that the matrix measures the current 
across one of the junctions instead of the number of electrons in the central region measured 
by N. 

In analogy with Eqs. (16) and (17) the current-current correlation function should have 
two matrices and for measuring current at the two times: and t. Although such 
a term includes the correlation between two different tunneling events, it does not include 
the self-correlation of a given tunneling event with itself. In our earlier paper-*^^ we gave 
a formal derivation of the self-correlation term as well as the correlation between different 
tunneling events. Here we put in this self-correlation term by hand. Suppose there is a 
current pulse across the right junction between and dt. The number current is {dt)~^ 
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during the time interval dt, and the number current squared is {dt) . The probabihty of 
this happening dt times the average number current to the right. Thus, it would seem that 
the self-correlation term for (/fl(t)/ii(0)) is just eI5{t). However, the current can be either 
positive or negative, and the self-correlation term is always positive. Thus, in computing 
the self-correlation term we must take the absolute value of the matrix elements of 
and used in computing the current. These new matrices may be written succintly as 
([N, v^])j j = |(v^)ij| and ([v^, N])jj = \{'v^)ij\- Thus, the current-current correlation 
function is {A, B — R or L) 

{lA{t)lBm = e2^(t)Tt{v^P(t)vV°'} 

+ e2^(-0lV{v^P(-0v V°^} (18) 

The first two lines of this equation contain the correlation between two different tunneling 
events, while the last line is the self-correlation term. 

As for the average current, the continuity equation identity, Eq. (13), can be used to 
show that any occurrence of {II — Ir) may be replaced by —e{dN/dt). 

-e|(Ar(t)iV(0)) = {{hit) - Ut))NiO)) (19) 

-ej^mt)iAm = {{hit) - m)hm. (20) 

In the zero frequency limit this implies that both the current-current correlation function 
and the current-density correlation functions are independent of where the current is mea- 
sured. The noise, Sab{<-^), is now defined as the fourier transform of the current-current 
correlation function, (/^(t)/B(0)), with its long time behavior subtracted off. 

/oo , . 

dte'^'UilAmsm-ilf). (21) 
-00 ^ ^ 

The continuity equation (Eq. (20)) and the relation ^'^^(a;) = S'ba(— <^) imply that cor- 
relation functions for all possible choices of A, B = R, L are equal in the zero frequency 
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limit, a; = 0. As discussed in the introduction we will always be taking the zero frequency 
limit because it is the easiest to study experimentally. 



E. Numerical Technique 



In order to compute the noise it is useful to symmetrize the matrix M. The key 
observation is that this system obeys detailed balance (see Eq. (8)): 

^i^jpr=^j^ipf- (22) 

A direct consequence of Eq. (22) is that M^jp^j^ = Mjipf\ allowing us to define a sym- 
metric matrix M, 



^ij ^ \Jpf/pT = V^ij^ir (23) 
In a similar manner we define v'^ as {p^j^ / pf^)^^"^ and pf^ as the square root of pf^ so 
that the expectation values may also be written in a symmetric form, e.g. 

= -e^(°'vV^- (24) 
The matrix M is written in terms of its eigenvalues, A < 0, and eigenvectors, xx, as 

Mij = J2Kxx)ii^x)r (25) 

A 

Using this same representation, the noise, -5^3(0;), is given by 

S^^iu) =2e2 y](p(o)v^x^)— 1— (x^vVO 

+2e2 y (pWv^x^)— ^ (x^vVO (26) 

In this paper we truncate the matrix M to include some finite number of states and then 
diagonalize it to determine its eigenvalues and eigenvectors. The eigenvector with zero 
eigenvalue is p'"^. Equations (24) and (26) are used to compute the current and noise. 
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III. LIMITS 



Before discussing the results of our numerical evaluation of the zero frequency noise, 
in this section we compute the noise in four limits which may be treated analytically. The 
cases discussed are: A. the zero voltage limit, where the noise is related to the conductance 
via the fluctuation dissipation theorem; B. the thermally activated conduction regime, 
where the voltage drop is large compared to the temperature but still small compared to 
the charging energy; C. the two state region, where the voltage is larger than the charging 
energy but small enough that only two charge states are allowed; and D. the large voltage 
bias limit. 

A. Zero voltage 

Even though we will focus on potential drops, eV, large compared to the thermal 
energy, fc^T, any technique for computing the noise must reproduce the fluctuation dissi- 
pation theorem^^'^^ at zero bias: 

S = 4kBTG{V = 0), (27) 

where G is the differential conductance, G = dl/dV. In the appendix we show how Eq. 
(27) follows from the definitions of the rates (Eq. (4)) and the expression for the noise (Eq. 
(26)). In Sec. IV. we also verify that our numerical algorithm produces the fluctuation 
dissipation theorem at zero bias. 

B. Thermally activated conduction 

When the voltage is large compared to the temperature, the thermal noise is small 
compared to the shot noise, which comes from the current flowing. Even with e|F| ^ ksT, 
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if the voltage is small compared to charging energy, the conduction is thermally activated 
because the voltage drop alone does not provide enough energy to overcome the charging 
energy. In this section we show that in this regime {Ec ^ e\V\ ksT) the noise is related 
to the current by the standard shot noise relation 



,5 = 2e|/|. 



(28) 



In this limit the steady state probability distribution function is strongly peaked 
about one value of n except for special values of Vp where two states are equally occupied. 
We shall assume that we are not at one of these special points and call the state which is 
most likely occupied the n = state. To describe the transport it is sufficient to keep only 
three states: the n = state and the two states immediately accessible from this state: 
n = ±1. At T = the matrix M thus becomes: 



M(T = 0) = 



-r_i^o 
r_i^o Ti^o 
-Fi^o- 



(29) 



The eigenvector of M(T = 0) with zero eigenvalue is (0,1,0), indicating that at zero 
temperature only p^_o is nonzero. For a finite temperature there is a small correction 
tfM = M — M(T = 0), which allows transitions to the states n = ±1. 



, -5r_i^o <^ro^-i 

SM = I 6T_,^o -STo^^ - STo^_, 6T,^o 
^To^i -(JTi^o, 



(30) 



There are analogous zero temperature and finite-temperature correction terms for and 



Because we have reduced M to a 3 x 3 matrix, equations (24) and (26) for the current 
and noise can easily be expanded to linear order in the ST^s: 



I 

e 

_S_ 
2? 



Ti^o 



-1^0 



1 + 



>0 



0^-1 



-1^0 



+ 



ri_ 



1^0 



►0 



-1^0 



(31) 
(32) 
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These equations have a simple physical interpretation. Conduction in this thermally ac- 
tivated regime follows from the occurrence of a rare transition from n = Oton = ±l, 
followed by a rapid decay back into the n = state. Thus, for example, there could be 
a thermal fluctuation causing an electron to tunnel from the right lead to the central re- 
gion (SFg^^). This thermal fluctuation is followed shortly by the electron tunneling off to 
the right lead (rf_^o) or to the left lead (rf_^o). If an electron tunnels to the right lead, 
then there is no contribution to the current, while if the electron tunnels to the left lead 
there is a positive contribution to the electrical current. Thus, only a fraction of the time, 
T^^q/Ti^o, does this thermally activated process contribute to the current. This explains 
the flrst term in Eq. (31): (rf_^g/ri^o)5r^_^^. The other terms in this equation have similar 
interpretations. 

The noise in this regime only contains the self- correlation term because we have inde- 
pendent thermally activated processes which are separated by long periods in time. Thus, 
the only difference between the current and the noise is one sign change resulting from 
some processes giving a negative contribution to the current and a positive self-correlation 
term for the noise. Dividing Eq. (32) by Eq. (31), it would seem that their ratio is not 2e; 
however, because the voltage is large compared to the temperature, two of the activated 
rates are always much larger than the other two. For / > 0, 5Tq_^-^^ and 5Tq^_^ are much 
larger than ST^^-^ and 5T^_^_^ by a factor exponential in e\V\/kBT, while for / < the 
situation is reversed. Therefore, the noise is indeed equal to 2e times absolute value of the 
current (Eq. (28)). 

C. Two state regime 

Eventually the voltage becomes larger than the charging energy, and there is current 
flow even at zero temperature. As discussed above, in the thermally activated regime only 
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one charge state is allowed at zero temperature. At the onset of non-thermally activated 
current flow, there are two allowed charge states at T = 0. It is a good approximation 
to keep only these two states for some range in temperature, ksT <^ Ec- Our two-state 
approximation is similar but not identical to the one used in Ref. 26. They compute p^"^ 
using only two states, and then allow transitions to a third state in computing the current 
(Eq. (24)). We do not allow transitions to a third state; however, in the regime of interest 
here, where only two states are energetically accessible, the two approaches give the same 
result for the current. 

Since only two states are involved, there are only two rates. In the case 7 > one 
rate, F^, increases the number of electrons by tunneling across the right junction, while 
the other rate, F^, decreases the number of electrons by tunneling across the left junction. 
These rates depend on the voltage drop across the sample and can be determined from 
Eqs. (3) and (4). In the next section we will also give the limiting form of these rates at 
zero temperature. Since M in this limit is a 2 x 2 matrix, 

^=(~X -r^V (33) 



^ F« -F 

it is simple to solve for the current and the noise 



I = e 



S = 2el 



(F^)2 + (F^)2 



(34) 

V ^ ' 



(35) 



(F« + F^)2 

We wish to emphasize here that although Eqs. (34) and (35) look simple, they contain 
a great deal of structure because the rates F^ and F^ depend on voltage. This will be 
shown explicitly in Sec. IV. when we plot this result along with the numerical results 
for finite temperatures. From Eqs. (34) and (35) we can see that / is largest when two 
rates are equal. Also, the noise is suppressed from its uncorrelated value, S — 2e|/|. This 
suppression is to be expected because having only two allowed states in the central region 
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introduces correlations between tunneling events. If the system is at one of the allowed 
states, say the state with more electrons, then another electron can not tunnel onto the 
central region until an electron tunnels off. This kind of correlation did not appear in the 
thermally activated regime because an electron tunnels off so quickly after tunneling on 
that it has no effect on subsequent tunneling events. 



D. High voltage limit 



We have now computed the noise for V = 0, Ec ^ e\V\ ^ ksT, and e\V\ ^ Eq- In 
this section we compute the noise for e\V\ 3> Ec-, where we expect the current- voltage 
characteristic to be almost linear. Even though we compute the noise in the asymptotic 
high voltage limit, in Sec. IV. we will see some of the approximations derived here work 
down to e\V\ ~ Ec- 

Because we are interested in the limit where the voltage is larger than the other energy 
scales, ksT and Ec, we set the temperature equal to zero. This means that 7 (e) in Eq. 
(3) takes on a very simple form: 7 (e) = e ^(e). If the sample is biased so that the electrical 
current flows from left to right, then two kinds of processes are important at T = 0: 
tunneling across the right junction to the central region {T^^^_^-^) and tunneling from the 
central region to the left lead {T^^^_-^). The rates for these processes are 

where the numbers and A^^ are given by 

iV. = -^-^ + i (38) 

C^I.C^ 1 

e e 2 
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These rates are illustrated in Fig. l.(b). The minimum allowed state, Ny^inj is the largest 
value of n for which (n — N^) < 0. Similarly, the maximum allowed state, N-max, is the 
smallest n for which {Nn — n) < 0. Assuming that Nl < and Nn > 0, this means that 
the values of n with are 

N^iri = -Int(-7Vi - 1) < n < Int(7V^ + 1) = Nmax- (40) 

Since In and 1^ are just the expectation values of T^^^_^-^ and T^_^^_-^ respectively, 
it is tempting to call Eqs. (36) and (37) with n — > (n) the current. This is not correct 
because T^^^_^ = {n — N^) / (RlC) is only valid for n > Nj^^i^ + 1- For n = N^^^ the rate 
i^ Similarly, T^^^^^ = {Nr - n)/{RRC) is only vahd for n < Nmax - 1- 

However, in the limit Rl <S Rr the rate for tunneling across the left junction, r^_^^_-,^, is 
much larger than the rate for tunneling across the right junction as illustrated in Fig. 1(b). 
This means that in) is close to Nmin- Because p^^ fa 0, it is then a good approximation 
to replace n by (n) in Eq. (37): 

^^e-^(Ar^-(n)). (41) 

In similar manner for R^ <^ Rj^ the probability of being in the state n = N^in is small 
and we can approximate the current by Eq. (36) with n —>■ (n): 

I ^ e^{{n) - N,). (42) 

Equations (41) and (42) show that the current measures the average number of electrons 
in the middle region. 

Up to this point we have not made any assumptions about the voltage being much 
larger than the charging energy. Now we assume that and Nj^ are integers. This can 
only be true for a limited set of voltages; however, in the large voltage limit the current 
and noise should insensitive to whether A^^ and A^^ are integers. We verify this numerically 
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in the Sec. IV. by showing that the results predicted here are asymptotically true in the 
high voltage limit. 

If Ni^ and are integers, then we have shown in our earlier paper-*^^ that the current 
and noise can be computed exactly. Here, we review the results in the context of this 
problem. For integer and Nji the number electrons in the middle region, N{t), satisfies 
a simple rate equation 

dNjt) ^ Nn - N{t) _ N{t) - 
dt RrC RlC ' ^ ' 

This allows one to compute both the steady state current 

and to relate the zero-frequency noise to the variance in the number of electrons in the 
middle region, var(n), 

S = 2el-— var(n) (45) 

Rd 



The variance of n is given by 



BZ'B;' (CV 



SO the noise in the high voltage limit is 



S = 2el , \ (47) 

Note the formal similarity of these results with those of the two state regime (Eqs. (34) 
and (35)). An important difference between these two sets of results is that R~j^ and Rr-x 
are constant, while and F^ depend on voltage. Thus, here S/2eI goes to a constant 
asymptotic value, while before S/2eI depended on voltage. In the next section we will 
show numerically that Eq. (47) is indeed valid in the high voltage limit, and Eq. (45) is a 
good approximation for much smaller voltages. 
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IV. NUMERICAL RESULTS 



In this section we present the results for computing the noise and current via Eqs. 
(24) and (26). There are six parameters in the model besides the voltage: the resistances, 
Rl and Rr, the capacitances, Cl and Cr, the temperature, T, and the offset voltage, Vp. 
Two of these parameters can be eliminated by rescaling the rates by i? = Rr + Rl and 
the energies by the charging energy, Eq = e^/2C. Thus, our parameter space consists 
of the resistance ratio, Rl/R, the capacitance ratio, Cl/C, the temperature divided by 
the charging energy, kBT/Ec, and eVp/Ec. In evaluating Eqs. (24) and (26) we will 
always truncate the matrix M because the higher charge states, \n\ ^ 1, are energetically 
forbidden, especially at low temperatures and voltages. In practice we have found it is 
sufficient to keep the 15 states surrounding n = 0. With this many states the steady state 
solution, is zero to within our numerical accuracy for the highest (n = 7) and lowest 
states (n = —7). Increasing the number of states does not change p^^. 



A. Temperature Dependence 



For most of this section we will consider temperatures far below the charging energy 
because this is where charging effects are most important. To start, however, we set the 
voltage equal to zero and verify that our numerical algorithm satisfies the fiuctuation 
dissipation theorem (Eq. (27)). In Fig. 2 we have plotted both the noise, 5, and the 
differential conductance, G = dl/dV, as a function of temperature for a pair junctions 
with Rl/R — 0.01 and Cl/C — 0.01. As expected in a coulomb blockade problem, 
the conductance and noise are suppressed for temperatures much less than the charging 
energy (^ 0.2£^(7, here). For temperatures greater than or of order the charging energy, 
the conductance rises linearly with temperature, indicating that the number of accessible 
states is increasing with temperature. In this same regime the noise approaches a constant 
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value as required by the fluctuation dissipation theorem. We have divided the two curves 
to show exphcitly that the fluctuation dissipation theorem is satisfled. 

B. Dependence on Vp 

In Fig. 3 the current and noise versus voltage curves are plotted for CVp = and 
CVp = 0.25. The temperature is much less than the charging energy, ksT = O.OlE^c- 
As expected for an asymmetric pair of junctions with = O.Oli? and = O.OIC, the 
curves exhibit the step-like structure of the coulomb staircase. The primary effect of the 
offset Vp is to shift the I-V and S-V curves. To reduce the parameter space in subsequent 
plots we will specialize to = 0, where the I-V and S-V characteristics are symmetric in 
the voltage. In this figure the fact that the noise and current curves are almost the same 
indicates that the noise is close to its uncorrelated value: S — 2el. Careful examination 
shows that the noise is actually less than 2el. The amount by which it is less than 2el 
changes as a function of voltage. Thus, in subsequent plots we will also look at the noise 
ratio, S/2eI, to bring out this structure. 

C. Dependence on the capacitance and resistance ratios 

Restricting ourselves to /c^T <^ Ec and Vp — 0, the parameter space becomes two 
dimensional: {Rl/Rr,Cl/Cr) or (Rl/R,Cl/C). It is well known that for asymmetric 
junctions with R^ <S Rr and <S Cr (or R^ » Rr and Cl ^ Cr) the coulomb staircase 
structure as in Fig. 3 is seen, while for asymmetric junctions with R^ ^ Rr and Cl ^ Cr 
(or Rl ^ Rr and C^ <^ Cr) no step-like structure is seen. There is only a suppression of 
the current for small voltages. 

These two limits can be understood as follows. The resistances, R^ and Rr, determine 
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the overall scale for the transition rates given in Eq. (4). If, for example, Rl <^ Rr^ then 
the tunneling rate across the left junction is much larger than the tunneling rate across 
the right junction. For a positive bias, > 0, electrons will tunnel off the middle region 
to the left lead. Thus, in the case Rl <S Rr-, the number of electrons in the central region 
will be close to the minimum allowed number, A^rnin of Eq. (40) . Using the approximate 
expression for the current when R^ <^ Rr, Eq. (41), there are steps in the I-V characteristic 
when (n) jumps, but no steps when (n) is approximately constant. 

This raises the question: when does (n) jump discretely with increasing voltage? While 
the resistances determine the overall scale of the rates, the capacitances determine the 
energetics. If Cl "C Cr, then the minimum allowed charge state at low temperatures is 
^min = -Int(CflF/e + 1/2) (see Eqs. (38) and (40)). On the other hand if Cl > Cr the 
minimum allowed charge state at low temperatures is N^^^^ = 0. Thus, in the first case 
the minimum allowed charge state changes with voltage while in the second case it does 
not. This means that for Cl <S Cr and Rl <i Rr the charge state will increase discretely 
as one increases the voltage, while for 3> Cr and Rl -C Rr the charge state will remain 
approximately constant. 

We have only discussed the two extreme limits when there are well defined steps in the 
I-V characteristic and when there are no steps in the I-V characteristic. Clearly, one can 
go continuously between these two cases by varying the resistance and capacitance ratios. 
In Fig. 4 we have computed the current and noise ratio, S/ el, as a function of voltage for 
Cl/C = 0.01 varying Rl/R, and in Fig. 5 we show the same quantities for Rl/R = 0.01 
varying Cl/C. For both ksT — O.lEc and /c^T = O.OlEc the I-V characteristic does 
indeed go between a step like structure for Rl <^ Rr and Cl <^Cr and an offset only for 
Rl <S Rr and Cl ^ Cr, or Rl » Rr and Cl <S Cr. In comparing Figs. 4 and 5 it is 
important to note that although the dimensionless ratio in the abscissa is the same in the 
two figures, the range of voltages is different. In Fig. 4 Cr is equal to 0.99C, while Cr 
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varies from 0.99C to 0.05C in Fig. 5. Thus, the curves at the top of Figs. 4 and 5 appear 
different, while they are in fact almost the same when plotted in terms of the unsealed 
voltage. 

Another important difference between Figs. 4 and 5 is the way one goes between the 
two extremes. When varying the resistances (Fig. 4), the steps become smoother, while 
in varying the capacitances (Fig. 5), the steps remain sharp, but the region between them 
gets a finite slope. In the first case when the resistances are changed, the probability 
distribution function becomes more spread out in n as becomes comparable to Rr. 
Thus, instead of having the distribution function peaked near A^rnin one obtains a much 
broader distribution function. Consequently, as new charge states become allowed, the 
distribution gradually shifts to lower n instead of jumping immediately to the new N^i^. 
In the second case when the capacitances are changed, the distribution remains fixed near 
the lower charge state, but the (n) in the current, / ~ {Nr — {n))/RRC, becomes less 
important relative to the term. In this case the jumps remain sharp, but become 
smaller in magnitude. 

Next we turn to the noise ratio S/2eI vs. current curves, which are the new contribu- 
tion of this paper. The first observation to make is that there is a rich variety of structure 
in the noise when rescaled by the current. This structure is not contained in the current. 
For example, one might say that the S/2eI curves in Fig. 5 are simply related to the 
derivative of the current; however, there is no simple way to take the derivative of the 
current curves in Fig. 4 and obtain the S/2eI curves in that figure. 

Using the structure in S/2eI, one can determine the parameters in the model more 
accurately than one could with the current alone. Alternately, one can use the parameters 
obtained from the current to perform a consistency check. We illustrate this point by 
considering three examples from Figs. 4 and 5. In Fig. 4(b), the top three cases, Rl/R = 
0.75, 0.9, and 0.99, have almost identical current- voltage characteristics, yet their noise 
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ratio curves are quite distinguishable. As a second example, the steps in the Rl/R = 0.5 
I-V curve are only slightly more pronounced than the steps in the Rl/R = 0.75 curve; 
however, the noise-ratio curves exhibit a change in curvature. The upper curve has a 
downward cusp, while the lower curve has an upward cusp. Finally, we note that the 
Rl/R — 0.1 curve in Fig. 4(b) and the CL/iCL + Cn) — 0.25 or 0.5 curves in Fig. 5(b) are 
quite similar, although their noise ratio curves are very different. 

Comparing Fig. 4(a) to 4(b) and Fig. 5(a) to 5(b), there is considerable sharpening in 
the noise ratio curves as one goes down in temperature. The spiky structure of the curves 
in Fig. 5(b) (/cbT = O.OlEc) is just beginning to be visible in Fig. 5(a) (ksT = O.lEc)- It 
is natural to ask whether even T = O.OlEc/kB is the zero temperature limit. For the higher 
voltages and less asymmetric junctions, one has indeed reached the zero temperature limit, 
i.e., the curves do not change appreciably as one goes down in temperature; however, the 
lower voltage regions of the more asymmetric junctions continue to become sharper as one 
goes down in temperature. To illustrate this we have plotted in Fig. 6 the region around 
the first step in the most asymmetric junction of Figs. 4 and 5 {Cl = O.OIC, Rl = O.Oli?) 
for four temperatures. The T — > curves were computed analytically using the two state 
approximation of Sec. III. C. For all of the curves shown, /, S, S/2eI, and dl/dV, there 
is significant enhancement of the structure for temperatures below one hundredth the 
charging energy. Both the current and the noise become sharper as one goes down in 
temperature, showing that the increased structure in the noise ratio is a combination of 
effects in the noise and current. 



D. Discussion of the noise ratio 



While the above shows that the noise provides new information about the transport in 
these systems, it does not give us a simple interpretation of what this information means. 
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In this section we use the analytic results of Sec. III. to understand our numerical results. 
Three limits are discussed: the thermally activated regime before the first step, the region 
just after the first step, where the two-state approximation is valid, and the large voltage 
limit, where the current and noise- voltage characteristics become linear. We also test an 
approximate formula for the noise which interpolates between the high and low voltage 
limits. This formula relates the zero frequency noise to the variance in the number of 
electrons in the middle region. All these cases except for the very high voltage limit are 
presented in Fig. 7. The asymptotic large voltage limit is illustrated in Figs. 4 and 5. 

The first region we consider is the thermally activated regime. In Sec. III. B. we 
showed that in this region S/2eI is unity because the current flows in one direction via 
rare uncorrelated events. In Figs. 4 and 5 we do indeed find that the noise ratio is unity 
in the low voltage regime below the first step in the I-V characteristic. This is shown 
more explicitly in Fig. 7, where we replot the solid curves of Figs. 4(b) and 5(b) along 
with several analytic approximations, shown as solid dots. In Fig. 7(a) the thermally 
activated regime extends up to CuV/e — 0.5 (region I). Because the capacitance ratios 
vary in Fig. 7(b), the extent of the thermally activated regime also varies. For Cr > Cl 
the thermally activated regime extends up to CuV/e = 1/2, while for Cl > Cr it extends 
up to CnV/e = (1/2){Cr/Cl)- Thus, the thermally activated regime for the top curve in 
this figure only extends up to CnV/e = 0.5(0.05/0.95) 0.026. We have not put in any 
solid dots for the thermally activated regime in Fig. 7(b) to avoid cluttering the graph. 
In both of these figures the noise ratio is not shown down to zero voltage because at any 
finite temperature {ksT = 0.01£?c here) thermal fiuctuations cause S/2eI to diverge as 
y ^ 0. 

The region just above the thermally activated regime in Fig. 7 is the two state limit. 
In this regime only two states have nonzero occupation probability. If the system is in one 
state, then the only allowed transition is to the other state. This introduces correlations 
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which in turn suppress the noise. For the processes in the thermally activated regime there 
were also two relevant rates: one slow rate due to thermal activation and a quick decay 
rate from the excited state. Because one of the rates was much longer than the other, the 
problem reduced to one with just the slow rate. Another way to say this is that the noise 
in the two state limit, S = 2e/((r^)2 + (r«)2)/(r^ + r^)^ (Eq. (35)), reduces to classical 
shot noise, S = 2el, when one of the rates is much larger than the other. 

In Figs. 7(a) and (b) one of the rates is much smaller than the other near the onset, 
and the noise ratio is close to unity. As the voltage is increased the smaller rate increases 
and the noise is suppressed. In this regime the maximum possible suppression is 1 /2 when 
pfl _ pi j^g^^ most of the curves in the two state regime (region II) do indeed 

go down to 1/2. The curves then return up as the rate which was originally smaller now 
becomes the larger rate. The exceptions are the Rl = O.Oli? and Cl = O.OIC junction, 
which has not yet reached the zero temperature limit (see Fig. 6), and the three upper 
curves, where the smaller rate does not reach the larger rate before the next step. In Fig. 
7(b) the same physics holds, but we have not illustrated the two state regime with solid 
dots because the place where the two state regime is valid varies with Cl/C. 

It is tempting to say that our two state approximation for the noise works even after 
the first step, especially for asymmetric junctions where the probability is strongly peaked 
about either the maximum or minimum allowed state. While the strictly two state ap- 
proximation predicts the qualitative features of the noise-ratio voltage curves for higher 
order steps, it is not quantitative. In particular we can see from Figs. 4 and 5 that the 
noise ratio does not go down to one half after the first step. Our two state approximation 
predicts that the noise goes down to 1/2 for the higher steps as well. Rather than try to 
extend the two state approximation by, for example, including three states, we now turn 
to the large voltage limit. 

In Sec. III. D. we argued that in the large voltage limit the noise ratio should be given 
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by S/2eI = (i?^ ^ + Rji'^)/{R^ ^ + R^^)'^ (Eq. (47)). This expression is remarkably similar 
to the expression we found in the two state limit for the noise (T^ji^ — * RL(ji))j however, it 
is important to keep in mind that the high voltage limit is definitely not a limit where the 
two state approximation is valid. Rather, formulas with this structure appear in a wide 
variety of contexts. ^'^'^^'^^ We regard them as a ubiquitous but not universal high voltage 
limit when only two rates are important. To check this asymptotic limit the dotted lines 
in Figs. 4 and 5 (not Fig. 7) are the ratio in Eq. (47). Clearly, it is approached at large 
biases. 

This still leaves the intermediate regime where more than two states are important 
and one has not reached the asymptotic limit of large voltages. From Eqs. (41) [R^ < Rr] 
and (42) [R^ < Ri] we see that at least for asymmetric junctions the current provides a 
measure of the mean number of electrons in the central region. In region III of Figs. 7(a) 
and (b) the solid dots over the numerical I-V characteristics show that these expressions 
are good approximations for the current even for the less asymmetric junctions. Since the 
current provides a measure of the mean number of electrons in the middle region, (n) , it is 
natural to try to express the current fluctuations in terms of the variance of n, var(n). In 
deriving the high voltage limit, we found that S ^ 2el - 4e\a.r{n) / {RC) (Eq. (45)). The 
solid dots in region III in Fig. 7(a) and all the solid dots in Fig. 7(b) are the noise ratio 
computed with this approximation using the numerically determined var(n). Although 
this is not as good an approximation as for the current, it works remarkably well down to 
lower voltages. The reductions in the noise can be regarded as an increase in the variance, 
which is to be expected near a step where there are large fluctuations. The places where 
the approximation works poorest, e.g. Rl/R = 0.5 and 0.75 in Fig. 7(a), are also the 
places where the approximation for the current work worst. This last approximation for 
the noise provides us with the simple intuitive picture that the noise measures the variance 
in n while the current measures the mean. 
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V. CONCLUSION 



In this paper we have computed the zero frequency current noise for the simple coufomb 
bfockade problem consisting of two tunnel junctions connected in series. The noise was 
computed both numerically and analytically using the same master equation as for the 
current. The numerical results are exact within the context of the master equation. The 
analytic results were shown to agree with the numerical results in four regimes: the zero 
voltage limit, the thermally activated regime, the region where only two charge states 
are accessible at low temperatures, and the large voltage limit. The low temperature 
limit in some cases is not obtained until the temperature is very much less than Ec, e.g. 
ksT = 10~^Ec. For the intermediate voltage regime we obtained a useful approximation 
to the noise which related the noise to the variance in the number of electrons in the 
region between the tunnel junctions. Thus, while the current measures the mean number 
of electrons in this region, the noise measures the variance. 

Both the numerical and analytic results showed that the noise contains information 
which is not contained in the current. Thus, by measuring the zero frequency noise we 
can determine the five parameters in the model more accurately than one could with the 
current alone. Alternately, one can use the noise to check the results obtained by measuring 
the current, showing possible deficiencies in the underlying rate equation used to describe 
transport in this system. 
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APPENDIX 



The proof of the fluctuation dissipation theorem in quantum mechanics is few hues 
long. On the other hand, the proof for most classical systems is much longer. In this 
appendix we prove the fluctuation dissipation theorem in two steps. We first simplify 
the expression for the noise in equilibrium and then simplify the expression for the linear 
response conductance. 

Simplify expression for the noise : Since we are interested in the zero frequency noise, 
we define a new matrix, SP, 

5Pij= r dt {Pij{t) - pf) . (48) 
J 

As can be seen by applying M to Eq. (48), 5Pij is related to the inverse of M: 

(M(JP),, = pf - 5ij. (49) 

The zero frequency noise may be expressed directly in terms SP. 

5^B = 2e2Tr{v^<JPvV°^} 

+ 2e2Tr{v^(5PvV°^} (50) 

In equilibrium the matrices and are proportional to one another. We thus define 
a matrix v which satisfies = -R^^v, = R^^, and 

^ ^N,M]. (51) 



Eq. (51) allows us to reduce Eq. (50) to two different kinds of traces: Tr{[N, [N, M]]p(°)} 
and Tr{[N, M]5P[N, M]p(°)}. These traces can be simplified using: (i) the trace of M 
times any vector is zero (Eq. (7)), (ii) Mp*^"^ is zero (Eq. (8)), and (iii) M5P is related to 
the identity matrix via Eq. (49). The end result is that the noise in equilibrium is 

^ = {R%%)2 ^{NMNp(°)}. (52) 
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Simplify expression for the linear response conductance : In equilibrium the probabil- 
ity function, has a thermal distribution, i.e. it is the exponential of — /3 times some 
energy. En. The simplest way to see this is to use detailed balance (Eq. (22)) 



n—l^n Pn e ^ " 

where the energy En is 



(53) 



2 2 
en 

= ^ + eV,n. (54) 

The fact that the distribution is thermal and that En is given by Eq. (54) allows one to 
express the derivative of p^°^ with respect to Vp in terms of the equilibrium probability: 



= -(3e (n - Tr{Np(°)}p(,°)) . (55) 



The strategy now is to express partial derivatives with respect to V in terms of partial 
derivatives with respect to Vp. At the end we use Eq. (55) to rewrite partials with respect 
to Vp in terms of equilibrium expectation values. Because we know that the currents II 
and are equal, we will only consider J^. The quantity which we are computing is the 
linear response conductance: 

^ = _eTV{^pW} - eTVlv^^}, (56) 
dV ^ dV'^ ^ ^ dV ^ ' 

where this and all other partial derivatives are evaluated at F = 0. The derivative of / 
with respect to at F = is zero because / at F = is zero. 

wr-^'^^wf^-^'^^'-k^^' (57) 

Because the F^'s are functions of CrV + CVp, derivatives with respect to V can be related 
to ones with respect to Vp. 

~dV ~ ~CdVp ^^^^ 
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Using Eqs. (56) - (58), one partial derivative with respect to V can be eliminated: 

In order to eliminate the other derivative with respect to V, we note that Mp^"' is zero 
for all voltages and hence 

As in Eq. (58), we can relate derivatives with respect to V to derivatives with respect to 



V,: 



(62) 



Using Eqs. (60) -(62), the linear response conductance is 

As our final simplification of dl/dV, we use Eq. (55) to covert the partial with respect to 
Vp to traces involving p'''^ 

§ = -"^^ I^{NMN.<»'} (64) 

Dividing Eq. (52) by Eq. (64) we obtain the the fluctuation dissipation theorem, Eq. (27). 
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FIGURE CAPTIONS 



Fig. 1 Schematic of experimental geometry and rates, (a) The capacitances and resistances of 
the left and right junctions are {Cl, Cr) and {Rl, Rr), respectively. The voltage in the 
left and right leads are fixed to V and 0, respectively. The voltage in the central region, 
Vm, fluctuates depending on the number of electrons in the central region, n. (b) At 
zero temperature the rates for tunneling are linear in n. For a positive bias, V > 0, 
only the rates for tunneling onto the central region from the right lead, r^^^_|_^, and 
the rates for tunneling off the central region to the left lead, T^^^_-^ are nonzero. The 
points where these rates vanish, A^^ and A^^, determine the maximum and minimum 
charge states, Nmax and iVrnin- 

Fig. 2 The equilibrium noise, S, and conductance, G = dl/dV, for two tunnel junctions 
connected in series. At temperatures far below the charging energy of the tunnel 
junctions, Ec = e^/2C, both the conductance (solid curve) and the zero-frequency 
noise (dotted curve) are suppressed by the coulomb blockade. At temperatures above 
the charging energy the conductance is roughly linear in temperature and the noise 
approaches a constant value. In accordance with the fluctuation dissipation theorem 
the ratio of the noise to the conductance remains fixed at 4/cbT (see dashed curve). 

Fig. 3 The nonequilibrium noise and current. At low temperatures {ksT = O.OlEc) for an 
asymmetric junction {R^/R = 0.01 and Cl/C = 0.01) both the average current, /, and 
the noise, S, show the characteristic step structure of the coulomb staircase. Indeed, 
S is roughly equal to the standard shot noise result of 2el. With an offset voltage, Vp, 
the current-voltage and noise-voltage characteristics are shifted, but other qualitative 
features of the curves remain the same. Henceforth we set = 0. 

Fig. 4 The dependence of the current and noise ratio {S/2eI) on the resistances of the 
junctions at (a) fc^T = O.lEc and (b) ksT = O.OlEc- We go from good steps 
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{Cl/C = 0.01,Rl/R = 0.01) to no steps {Cl/C = 0.01, Rl/R = 0.99) by vary- 
ing the resistance ratio, Rl/Rr, keeping the capacitances fixed at Cl/C = 0.01 and 
Cr/C = 0.99. When the difference between S and / iUustrated in Fig. 3 is plotted as 
the noise ratio, S/2eI, a rich variety of structure is revealed. This structure becomes 
more pronounced as one goes down in temperature. The noise ratio shows structure 
even when the current shows little structure, illustrating how the noise can provide new 
information about the parameters in the model. The dotted lines are the asymptotic 
value of the noise ratio at large voltages (Eq. (47)). Since the noise remains finite as 
the current goes to zero, the noise ratio diverges as F — > 0. We have thus only plotted 
it for CrV/c > 0.2. (This is the reason S/2eI > 1 at low voltages in (a).) 

Fig. 5 The current and noise ratio as a function of the capacitance ratio at (a) ksT = O.lEc 
and (b) ksT = O.OlEc- {Rl/R = 0.01) Varying the capacitance instead of the re- 
sistance has a different effect. All the noise ratio curves have a similar structure, 
indicating that the noise is more sensitive to the the resistance ratio than capacitance 
ratio. This should not be too surprising because the asymptotic limit of noise ratio 
(dotted lines) is determined by the resistances of the junctions (see Eq. (47)). Al- 
though the abscissas of Fig. 4 and this figure are the same, the voltage ranges are 
different because the capacitances vary here. 

Fig. 6 The temperature dependence for the most asymmetric junctions shown in Figs. 4 and 
5 {Rl/R = 0.01, Cl/C = 0.01). The current (/), the noise (5), the noise ratio {S/2eI), 
and the differential conductance {G = dl/dV) all show significant temperature depen- 
dence below ksT = 0.01£^c- The spiky structure in the noise ratio is due to a combined 
effect of / and S having sharper structure as one goes to lower temperatures. The zero 
temperature curves are obtained from the analytic expressions for the noise in the two 
state limit (Eqs. (34) and (35)). 

Fig. 7 Approximations for the current and the noise. The solid lines are the numerical results 
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of Figs. 4(b) and 5(b), and the dots are the approximations, (a) In the thermally 
activated regime labeled by I the noise ratio is unity. In the two state regime labeled 
by II the current and noise are given by Eqs. (34) and (35). Both this and the thermally 
activated regime are exact at low temperatures. In the higher voltage regime, region 
III, we have used Eqs. (41) [Rl < Rr] -(42) [Rr < Rl] and Eq. (45) which relate 
the current and noise to the mean and variance in the number of electrons in the 
region between the two tunnel junctions. Although these approximations are not exact, 
they clearly get the important qualitative features of the current and noise, (b) With 
this abscissa the boundaries of the thermally activated and two-state regimes vary as 
one changes the capacitances. Thus, although the thermally activated and two state 
regimes are also exact here, we have only shown the large voltage approximation to 
avoid cluttering the graph. 
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